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"by  J.  J.  Stoker  and  A.  S.  Peters 


1*  Introduction  and  summary* 

The  purpose'  of  this  report  is  to  develop  the  mathematical 
theory  for  tho  motion  of  a ship,  to  bo  treated  as  a freely  floating 
rigid  body  under  the  action  of  given  external  forces  (a  propeller 
thrust,  for  example),  under  tho  most  general  conditions  compatible 
with  a linoar  thoory  and  the  assumption  of  an  infinite  ocoan# 

This  of  course  requiros  the  amplitudo  cf  tho  surface  waves  to  be 
small  and,  in  gonoral,  that  tho  motion  of  tho  water  should  bo 
small  oscillations  near  its  rest  position  of  equilibrium*  it  also 
requires  the  ship  to  havo  tho  shape  cf  a thin  disk  so  that  it  may 
have  a translatory  motion  with  finite  velocity  and  still  create 
only  small  disturbances  in  tho  water.  In  addition,  tho  motion  of 
tho  ship  itsolf  must  bo  assumed  to  consist  of  small  oscillations 
relativo  to  a uniform  translation.  Within  theso  limitations, 
however,  tho  thoory  to  bo  prosontod  is  quite  gonoral  in  tho  sonso 
that  no  arbitrary  assumptions  about  tho  interaction  of  tho  ship 
with  tho  wator  oro  made,  nor  about  the  character  of  tho  coupling 
betwoon  tho  difforont  degroos  cf  froodom  of  tho  ship,  nor  about 
tho  wavos  prosent  on  tho  surfaco  of  tho  soa;  tho  combined  system 
of  ship  and  sea  is  treated  by  using  tho  basic  mathematical  thoory 
of  tho  hydrodynamics  of  a non-turbulont  perfoct  fluid. 


For  example,  .the  theory  presented  here  would  make  it  possible  to 
determine  the  motion  of  a ship  under  given  forces  which  is  started 
with  arbitrary  initial  conditions  in  a sea  subjected  to  given  sur- 
face pressures  and  initial  conditions,  or  in  a sea  covered  with 
waves  of  prescribed  character  coming  from  infinity* 

It  is  of  course  well  known  that  such  a linear  theory  for 
the  non-turbulent  motion  of  a perfect  fluid,  complicated  though 
it  is,  still  does  not  contain  all  of  the  important  elomonts  needed 
for  a thoroughgoing  discussion  of  the  practical  problems  involved. 
For  example,  it  ignoros  tho  boundary-layer  effects,  turbulent 
effects,  the  existence  in  general  of  a wake,  and  other  important 
effects  of  a non  linear  character*  Good  discussions  of  theso 
matters  can  bo  found  in  papers  of  Lunde  and  Wigloy  [6]  and  Jfcvulook 
[3]«  Nevertheless,  it  sooms  dear  that  an  approach  to  tho  problem 
of  predicting  mathematically  tho  motion  of  ships  in  a seaway  under 
quite  general  conditions  is  a worth-whilo  enterprise,  and  that 
a start  should  bo  made  with  tho  problom  oven  though  it  is  rocog- 
nizod  at  tho  outset  that  all  of  the  important  physical  factors 
can  not  bo  taken  into  account.  In  fact,  tho  theory  presented  here 
loads  at  onoo  to  a number  of  important  qualitative  statements 
without  tho  necessity  of  produo ing  actual  solutions  - for  example, 
wo  shall  see  that  certain  resonant  frequencies  appear  quite 
naturally,  and  in  addition  that  thoy  con  bo  oaloulatod  solely 
with  roferonco  to  the  mass  distribution  and  tho  given  shapo  of 
tho  hull  of  the  ship*  Intorosting  observations  about  tho  character 
of  tho  coupling  between  tho  various  dogroos  of  freedom,  andebout 
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the  nature  of  the  interaction  between  the  ship  and  the  water,  are 
also  obtained  simply  by  examining  the  equations  Which  the  theory 
yields. 

In  order  to  describe  the  theory  arid  results  to  be  worked  out 
in  later  sections  of  this  report,  it  is  necessary  to  introduce  our 
notation  and  to  go  somewhat  into  dotr.ils.  In  Pig.  1.1  the  dispo- 
sition of  two  of  the  coordinate  systems  used  is  indicated.  The 


Pig . 1*1. 

Fixed  and  Moving  Coordinate  Systems 


system  (X,Y, Z)  is  a system  fixed  in  space  with  the  X,Z-plane  in 
the  undisturbed  free  surface  of  the  water  and  the  Y-axis  vortic  ally 
upward.  * 

ft 

This  choice  of  axis  is  not  the  conventional  one;  the  Z-axis  is 
usually  chosen  as  the  vortical  axis.  It  wa3  made  because  the 
authors  aro  accustomed  tc  working  with  a variety  of  different 
water  wave  problems;  and  the  choice  made  here  seemed  to  them  to 
bo  reasonable  from  a general  point  of  view  because  of  the  largo 
number  of  existing  two-dimensional  problems  of  interest  in  which 
one  then  naturally  chooses  the  y-axis  as  vertical  axis,  coupled 
with  the  fact  that  the  use  of  the  symbol  z as  a complex  variable 
is  nearly  universal. 


A moving  system  of  ooordinatos  (x,y,z)  is  introducod|  in 
this  system  the  x,z»plano  is  assumed  to  coincide  always  with  the 
X,Z-plone,  and  its  y-axis  is  assumed  to  contain  tho  center  of 
gravity  (abbroviatod  to  c.g.  in  tho  following)  of  tho  ship,  Tho 
courso  of  tho  ship  is  fixed  by  tho  motion  of  the  origin  of  tho 
moving  system!  it  is  then  convenient  to  introduce  the  speed  s(t) 
of  the  ship  in  its  course;  the  speed  s(t)  is  simply  tho  magnitude 
of  tho  vector  representing  tho  instantaneous  volocity  of  this 
point.  At  the  samo  time  wo  introduce  the  angular  spood  ui  ( t ) of 
tho  moving  systom  relative  to  tho  fixod  system*  ono  quantity  fixes 
this  rotation  because  tho  vortical  axes  remain  always  parallel. 

The  angle  a(t)  indicated  in  Pig, 1.1  is  defined  by 

t 

(1,1)  a(t)  = j u(t)dt  . 

0 

In  order  to  doal  with  tho  rigid  body  motion  of  tho  ship  it  is 
convenient,  as  always,  to  introduce  a systom  of  coordinates  fixod 
in  tho  body.  Such  a systom  (x*,y*,z»  ) is  indicatod  in  Pig.  1.2 


/ 

it?,' 

Fig,  1,2 

Tho  Moving  Coordinate  Svstom 


Tho  *' ,y*-plnno  Is  assumed  to  be  in  the  f oi’o-ond-af t plane  of 
symmetry  cf  the  ship*s  hull,  and  the  y»-nxis  is  assumed  to  contain 
tho  c.g,  of  tho  ship.  The  moving  systom  (x*,y',z*)  is  assumed  to 
coincido  with  tho  (x,y,z)  system  when  tho  ship  and  tho  water  are 
at  rost  in  thoir  equilibrium  positions.  The  c.g.  of  tho  ship  will 
thus  coincide  with  tho  origin  of  tho  (x',y*,z*)  system  only  in 
case  it  is  at  the  level  of  tho  equilibrium  water  line  on  the  shipj 
we  therefore  introduce  tho  constant  y*  as  tho  coordinate  cf  the 
c,g,  in  tho  primod  coordinate  systom  at  such  an  instant. 

^ of  tho  wator  is  assumed  to  be  given  by  a volocity 

potential*  f(X,Y,Zjt)$  it  in  turn  is  thoroforo  to  bo  determined  as 
a solution  of  Laplace* s equation  satisfying  appropriate  boundary 
conditions  at  tho  freo  surface  of  tho  wator,  on  tho  hull  of  tho 
ship,  at  infinity,  and  also  initial  conditions  at  tho  tiroo  t * 0. 
The  boundary  conditions  on  tho  hull  of  the  ship  clearly  will  depend 
on  the  motion  of  tho  ship,  which  in  its  turn  is  fixed,  through  tho 
differential  equations  for  the  motion  of  a rigid  body  with  six 
degrees  of  freedom,  by  tho  forces  acting  on  it  - including  tho 
prossuro  of  tho  water  - and  its  position  and  velocity  at  the  time 
t = 0.  As  wo  have  alroady  stated,  wo  make  no  further  restrictive 
assumptions  except  those  neodod  to  linearize  the  problem. 

Boforo  discussing  tho  linearization  wo  interpolate  a brief 
discussion  of  the  relation  of  the  presont  work  to  that  cf  cthor 
writers  who  have  discussed  tho  problem  of  ship  motions  by  means 

’^n-^soouaf?iiw.that  W°  d0al  Klth  0,1  ^rotational  motion  of  a 
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of  the  linear  theory  of  irrotational  waves.  Tho  subject  has  a 
lengthy  history,  beginning  with  Mioholl  [ 8 ] in  1898,  and  continu- 
ing over  a long  period  of  yoars  in  a sequence  of  notable  papers  by 
HaVolock,  boginning  in  1909.  This  work  is,  of  courso,  included  as 
^ a special  case  in  what  is  presented  hero.  Extensive  and  up-to-date 

bibliographies  can  bo  found  in  tho  papers  of  Weinblum  [10]  and 
Lundo[7].  Most  of  this  work  considers  tho  ship  to  bo  hold  fixod 
in  spaco  while  the  water  streams  past}  tho  question  of  intorost 
is  then  the  calculation  of  the  wave  resistance  in  it3  dependence 
on  tho  form  ef  the  ship.  Of  particular  intorost  to  us  here  are 
papors  of  Krylov  [5],  Weinblum  and  St.  Denis  [9],  and  Haskind.  [1], 
all  of  whom  deal  with  loss  restricted  typos  of  motion.  Krylav 
t seeks  tho  motion  of  the  ship  on  tho  assumption  that  the  prossuro 

on  its  hull  is  fixed  by  tho  prescribed  motion  of  tho  water  without 
roforonco  to  the  back  effect  on  tho  motion  of  tho  watc-r  induced  by 
the  motion  of  the  ship.  Weinblum  and  3t.  Denis  omplcy  a combined 
theoretical  and  empiric al  approach  to  tho  problem  which  inv<~lvos 
writing  down  equations  of  motion  of  tho  ship  with  coefficients 
which  should  bo  in  part  determined  by  model  experiments;  it  is 
assumed  in  addition  that  there  is  no  coupling  between  tho  different 
degroos  of  freedom  involvod  in  tho  general  motion  of  the  ship* 
Haskind  attacks  tho  problom  in  tho  some  dogroo  of  generality, 

. under  tho  snrno  general  assumptions,  as  tho  authors;  in  tho 

end,  however,  Haskind  derives  his  theory  completely  only  in  a 
certain  special  case.  Haskind1 s theory  is  also  not  the  sano  ns 


tho  thoory  presented  hore,  and  this  is  caused  by  a fundamental 
difference  In  tho  procedure  used  to  derive  tho  linear  thoory  from 
tho  underlying,  basically  nonlinear,  theory*  Haskind  dovolops  His 
thoory,  in  the  time -honored  way,  by  assuming  that  ho  knows  a jprl-or*! 
the  relative  orders  of  magnitude  of  tho  various  quantities  inwol.ve>d 
Applied  mathematicians  are  not  often  deceived  in  following  such  a 
procedure,  but  the  present  case  is  exceptional  both  because  of 
its  complexity  and  because  of  tho  fact  that  it  is  essential  to 
consider  terms  which  are  not  all  of  tho  same  order#  The  authors 
also  tried,  to  attack  the  problem  (without  being  aware  at  tho  feirao 
of  the  existence  of  Haskind' s work)  in  this  sarao  way,  but  invari- 
ably arrived  at  formulations  which  soemod  to  bo  inconsistent* 
Consequently  they  felt  it  necessary  to  proceed  by  a formal  devoLop- 
ment  with  rcspoct  to  a small  parameter  (essentially  tho  breadtli- 
longth  ratio  of  the  ship)}  in  doing  so  ovory  quantity  was  dove  looped 
systematic  ally  in  a formal  serios  (for  a similar  typo  of  discussion 
see  P*  John  [ijJ)*  In  this  way  a cerroct  theory  should  bo  obtained , 
assuming  tho  convergence  of  tho  serios  - and  tho  authors  see  no 
reason  to  doubt  that  tho  sorios  would  convorgo  for  sufficiently 
small  values  nf  the  parameter.  Asi'-’o  from  tho  relative  safety 
cf  such  a method  - purchased,  it  is  true,  at  tho  price  of  makixig 
rather  bulky  calculations  - it  has  an  additional  advantage,  i*e>,„ 
it  makes  possible  a consistent  procedure  for  determining  ray 
dosired  higher  order  corrections#  It  is  not  easy  to  compare 
Haskind* s theory  in  dotail  with  tho  thoory  prosontod  hero#  Howovor* 


it  con  bo  statod  that  cortain  terms,  oallod  damping  terms  by 
Haskind  and  oonsidorod  to  be  of  Importance  by  him,  arc  terms  that 
would,  in  the  theory  presented  hero,  bo  of  higher  order  than  any 
of  the so  rotained  by  the  authors}  consoquontly  tho  authors  fool 
that  conclusions  drawn  from  such  terms  may  well  bo  illusory  unloss 
somo  ovidon.ee  is  presented  which  shows  thoso  torms  to  bo  tho  most 
important  among  tho  very  largo  number  of  different  torms  of  that 
order  which  would  occur  in  a formal  dovolopraont.  A more  prociso 
statement  cn  this  point  will  be  made  later* 

Tho  procedure  followed  hero  begins  by  writing  the  equation 
of  tho  ship*s  hull  relative  to  tho  coordinate  system  fixed  in  tho 
ship  in  the  form 

d*2)  z*  = + ph(x*,y>)  , z»  ^ 0 , 

with  p a small  dimensionless  parnmotor*  This  is  tho  parameter 
refor  od  to  above  with  respect  to  which  all  quantities  will  bo 
developed*  In  particular,  tho  velocity  potential  4(X,Y,Z$tjp)  s 
/(*»y#**t|P)  is  assumod  to  possess  tho  development 

(1.3)  /*(x,y,zjt}p)  = p/^U.yjZjt)  + p2/rf£(x,y,z}t)  + •••••  , 

Tho  freo  surface  olovation  i\(x,zjt}p)  and  tho  speed  s(t}p)  and 
angular  velocity  u(t}P)  (cf.  1.1)  are  assumed  to  have  the 
developments 


(1.4) 

(1.5) 

(1.6) 


• • • 


I 


n(x,zjtjp)  * + p2n2(x,zjt)  + 

s(t$p)  * SQ  { t ) + ( t ) + ••• 

cj(t|p)  « ^c(t)  + pwL(t)  + ••• 

/■ 

Finally,  the  vertical  displacement  yc(t)  of  the  center  of  gravity 
and  the  angular  displacements  6^,  ©2,  of  th®  ship  with  respect 
to  the  x,y,  and  z axes  respectively  are  assumed  given  by 

(1.7)  ©x(t|p)  a p en(t)  + P2  ©12(t)  ♦ •••,  i a 1,2,3, 

(1.8)  y0(tjp)  * y£  * Py^t)  * P2y2(t)  + •••  • 

These  relations  imply  that  the  velocity  of  the  water  and  the 
elevation  of  its  free  surface  are  small  of  the  same  order  as  the 
"slenderness  parameter"  p of  the  ship#  On  the  other  hand,  the 
speed  s(t)  of  the  ship  is  assumed  to  bo  of  zoro  order*  The  other 
quantities  fixing  the  motion  of  tho  ship  are  assumed  to  bo  of . 
first  order,  oxcapt  for  u(t),  but  it  turns  out  in  the  end  that 
w0(t)  vanishes  so  that  u is  also  of  first  order*  Tho  quantity 
y£  in  (1*8)  was  dofined  in  connootion  with  tho  description  of 
Flg.1.2 j it  is  to  be  noted  that  wo  hnvo  chosen  to  oxpross  all 
quantities  with  rospoct  to  the  moving  coordinate  system  (x,y,z) 
Indicated  in  that  figuro*  Tho  formulas  for  changes  of  coordinates 
must  be  used,  and  thoy  also  are  to  bo  dovolopod  in  powers  of  pj 
for  example,  tho  equation  of  tho  hull  relative  to  the  (x,y,z) 
coordinate  system  is  found  to  be 
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z + p©21x  - P©JLl(y-y^)  - f*h(x,y)  + • * o 


after  developing  and  rejecting  seoond  and  highor  ordor  terms  in  p# 
In  merino  engineering  there  is  an  accepted  terminology  for 
describing  the  motion  of  a ship#  wo  wish  to  put  it  into  relation 
with  tho  notation  just  introduced#  The  angular  displacements  are 
named  as  follows l ©^  is  the  rolling,  fig+a  is  tho  yawing , and 
©j  is  the  pitching  oscillation#  The  quantity  s^(t)  in  (1 #5)  is 
called  the  surge  (i#o»,  it  is  tho  small  fore-and-aft  motion  rela- 
tive to  the  f ini to  spood  sQ(t)  of  tho  ship)#  while  y fixes  tho 
heave#  In  addition#  there  is  the  sido-wise  displacement  (in  first 
order  it  might  bo  denoted  by  pz^t))  roforred  to  as  the  sway  I thi3 
quantity,  in  lowest  ordor,  can  be  calculated  in  terms  of  sQ(t)  and 
the  anglo  a defined  by  (1,1)  in  terms  of  u>(t)  as  follows: 


(1.9) 


Mx(t) 


soa 


P«0  J ^ ( t )dt , 
0 


since  i.*o(t)  turns  out  to  vanish#  In  ono  of  the  problems  of  most 
practical  interest#  i#o#  the  problem  of  a Ehlp  that  has  been  moving 
for  a long  time  (so  that  all  transients  hr.vo  disappeared)  under 
a constant  propeller  thrust  (considered  to  bo  simply  a fcrco  of 
constant  magnitude  parallel  to  tho  keel  of  the  ship)  into  a soaway 
consisting  of  a given  system  of  simple  haxvnonic  progressing  wc.vos 
of  given  frequency#  ono  expects  that  the  displacement  components 
would  in  general  bo  the  sum  of  two  terms#  one  independent  of  tho 
time  and  representing  tho  displacomonts  that  would  arise  from 


f-.; 


-11- 


motion  with  uniform  volocity  through  a calm  sea,  tho  other  a term 
simple  harmonic  in  the  time  that  has  its  origin  in tho  foroos 
arising  from  tho  waves  coming  from  infinity.  On  account  of  the 
symmetry  of  tho  hull  only  two  displacements  of  tho  first  category 
would  differ  from  zeros  one  in  the  vortical  displacement,  i.e. 
tho  heave,  tho  other  in  the  pitching  angle,  i.o.  tho  angle  0^. 

Tho  lattor  two  displacements  apparently  are  referred  to  as  tho 
trim  of  tho  ship.  In  all,  then  there  would  bo  in  this  case  nine 
quantities  to  bo  flxod  as  far  as  tho  motion  of  tho  ship  is*  con- 
cerned: tho  amplitudes  of  tho  oscillations  in  each  of  tho  six 
dogroos  of  freedom,  tho  speed  sQ,  and  tho  two  quantities  deter- 
mining the  trim. 

Wo  procood  to  givo  a summary  of  the  theory  obtained  when  tho 
series  (1.2)  to  (1.8)  are  inserted  in  all  of  tho  equations  fixing 
the  motion  of  thG  systom,  which  includes  both  tho  differential 
equations  and  tho  boundary  conditions,  and  any  functions  involving 
P are  in  turn  dovolopod  in  powers  of  p.  For  oxamplo,  one  needs 
to  evaluate  on  the  free  surface  y * vj  in  ordor  to  oxpros3  tho 
boundary  conditions  there j one  calculatos.it  as  follows  (using 
(1.3)  and  (1.4)* 

* * * 

(1.10)  /x(x,^,zjtip)  * p[^lx(x,0,zjt)  + (x,0,zjt)  ♦ •••] 

• • * • • f 

+ P2^2x  + ^2xy  + **•!  * ••• 

■ P/l3C(x*0,zjt)  + P^Ui/ixyUfOjZjt)  + /^(XjOjZjt))  + •••  , 


,:*»**■ 
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Wo  obsorve  tho  Important  fact  - to  which  roforonce  will  be  made 
later  - that  the  coefficients  of  tho  powers  of  p are  ovaluatod 
at  y * 0,  l.o.  at  tho  undisturbed  oquilibrlum  position  of  tho 
free  surface  of  the  water#  Tho  oncl  rGsult  of  such  calculations, 
carried  out  in  such  a way  as  to  includo  all  torms  of  first  order 
in  P is  as  follows*  The  differential  equation  for  is,  of 
courso,  tho  Laplace  equation* 


(1.11) 


d.  + d.  + d-,  * 0 

^lxx  ^lyy  Flzz 


in  tho  domain  y < 0,  i«o#  the  lower  half-space,  excluding  tho 
plane  area  A of  the  x,y-plono  which  is  tho  orthogonal  projection 
of  tho  hull,  in  its  equilibrium  position,  on  tho  x,y-plnno#  The 
boundary  conditions  on  d}  ai>0 


(1.12) 


f *u  88  so(hx"°2i)  * (wi^2i)x  " on  K 

I 

i/lz  = “so(hx+®21)  + H+®21)x  " ®ll(y"ye)»  on  A- 


in which  A.  and  A rofor  to  the  two  sides  z - 0 . and  z * 0 of 

T • T • 

the  plane  disk* A.  Tho  boundary  conditions  on  tho  free  surface  are 


(1.13) 


8,lX  + aJlx  - *lt  * 0 

hv  • “o’U*  * ’'it  * 0 


at  y « 0, 


The  first  of  these  results  from  the  condition  that  the  press  re 
vanishes  on  the  free  surface,  the  second  arises  from  the  kinematic 
free  surface  condition.  If  soic»^,$2i#  iUlc*  ^11  were  known  functions 
of  t,  these  boundary  conditions  in  conjunction  with  (1.11)  and 
appropriate  initial  conditions  would  serve  to  determine  the  function 
and  n^  uniquely;  i.e.  the  velocity  potential  and  the  free  surface 
elevation  would  be  known.  In  any  case,  the  function  ^ - which  we 
repeat,  fixes  the  lowest  order  term  in  the  development  of  the 
velocity  potential  ^ - could  be  in  principle  determined,  because 
of  the  linearity  of  the  problem,  as  a linear  combination  of  harmonic 

functions  \|f^  having  s^w^  + ®21*®21  and  ®11  as  time“dePendont 
coefficients: 

(1.14)  ^(x.y^jt)  = s^U^y,*)  + (wj  + ®2i)*2*x'y,z^+*21%*x,y'z* 

+ ®nVx,y,z)  + 

The  harmonic  function  would  be  determined  through  initial  con- 
ditions and  tho  condition  fixing  the  wave  train  coming  from  od  - 
that  is,  it  contains  the  part  of  the  solution  arising  from  the 
nen-hdmogeneous  conditions  in  the  problem* 

Before  continuing  to  describe  the  relations  which  determine 
the  time -dependent  coefficients  in  (l.lq.)  as  well  as  tne  other 
unknown  functicns  of  the  time  which  fix  the  motion  of  the  ship,  we 
digress  at  this  point  in  order  to  discuss  some  conclusions  arising 
from  our  developments  and  concerning  two  questions  which  recur 
again  and  again  in  the  literature.  These  issues  center  around  tho 
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question:  what  is  tho  correct  manner  of  satisfying  the  boundary 
conditions  on  the  curved  hull  of  the  ship?  Michell  employed  tho 
condition  (1.12),  naturally  with  «n  * ©21  ■ ^ ■ 0,  on  the  basis 

of  the  physical  argument  that  sohx  represents  tho  compcnont  of  tho 
velocity  of  tho  water  normal  to  tho  hull  and  since  tho  hull  is 
slender,  a good  approximation  would  rc.  ult  by  using  as  boundary 
condition  the  Jump  condition  furnished  by  (1.12).  Havelock  and 
others  have  usually  followed  tho  same  practice.  However,  enu  finds 
constant  criticism  of  the  rosulting  theory  in  the  litoraturo 
(particularly  tho  engineering  litoraturo)  because  of  the  fact  that 
tho  boundary  condition  is  not  satisfied  at  the  actual  position  of 
tho  ship's  hull,  and  various  proposals  have  boon  mado  to  improve 
tho  approximation.  The  authors  fool  that  this  criticism  is  beside 
tho  point,  since  tho  condition  (1.12)  is  simply  tho  consequence  of 
a roasonnble  linearization  of  tho  problem.  To  take  account  of  tho 
boundary  condition  at  th<-  actual  position  of  tho  hull  would,  of 
course,  bo  more  accurate--  but  thon,  it  would  be  nocossary  to  deal 
with  tho  full  nonlinear  problem  and  mako  sure  that  all  of  tho 
essential  correction  terms  of  a given  order  wore  obtainod#  In 
particular,  it  would  be  necessary  to  examine  tho  higher  order 
terms  in  tho  conditions  at  tho  froo  surface  - after  all,  tho  con- 
ditions (1.13),  which  aro  also  usod  by  Michell  and  Havelock,  are 
satisfiod  at  y * 0 nnd  not  on  tho  actual  displaced  position  of  tho 
froo  aurfaoo.  Ono  way  to  obtain  a roar©  accurate  theory  would  bo, 
of  courso,  to  carry  out  tho  perturbation  scheme  outlined  hero  to 
higher  order  terms. 


v 
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Still  another  point  has  come  up  for  froquont  discussion  (cf#, 

for  oxomplo,  Lunde  and  Wigley  £63 ) with  roforenco  to  tho  boundary 

condition  on  tho  hull#  It  is  fairly  oommon  in  tho  literature  to 

refer  to  ships  of  Micholl‘s  typo,  by  which  is  meant  ships  which 

are  slender  not  only  in  tho  f oro— and— nf t direction,  but  which  are 

also  slondor  in  the  cross-sections  at  right  angles  to  thi3  tiiroc- 

tion  (cf#  Pig#l#3)  so  that  h^.,  in  our  notation*  is  small#  Thus 

%> 

ships  with  a rather  broad  bottom  (cf#  Fig#l#3),  or,  as  it  is  also- 
put,  with  a full  mid-section,  are  often  considered  as  ships  to 


bocomo  rather  largo  near  tho  kool  of  tho  ship  for  certain  typos 
of  cross-sections,  but  nevertheless  tho  linearization  carried  out 
above  should  remain  valid  sinoo  all  that  is  noedod  is  that  tho  ship 
should  not  oroato  too  largo  disturbances,  and  this  condition  is 
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guaranteed  by  taking  a long,  slondor  ship.  (It  might  also  be  noted 
that  h occurs  in  our  theory  only  under  integral  signs.)  In  fact, 

w 

there  c.ro  experimental  results  (cf.  Havelock  [3]  ) which  indicato 

that  the-  theory  is  just  as  accurate  for  ships  with  a full  mid* 

section  as  it  is  for  ships  of  Miclicll's  typo. 

After  this  digression  we  return  once  more  to  the  description 

of  the  equations  which  dotormino  the  motion  of  tho  ship,  and  which 

arise  from  developing  the  equations  of  motion  with  respect  to  0 

2 

and  retaining  -nly  tho  terms  of  order  p and  p . (Wo  observe  again 

that  it  is  necessary  to  consider  terms  of  both  orders. ) In  doing 

so  tho  mass  M of  tho  ship  is  givon  by  M = M^p,  with  a constant, 

since  we  assume  tho  averago  density  of  tho  ship  to  bo  f ini to  and 

its  volume  is  of  course  of  order  p.  Tho  moments  of  inertia  aro 

also  of  ordor  p.  Tho  propollor  thrust  is  assumed  to  bo  a forco  cf 

magnitudo  T acting  in  tho  negative  x* -direction  and  in  tho 

x* ,y' -piano  at  a point  whose  y* -coordinate  is  "JLi  tho  thrust  T is 
2 

cf  ordor  p , since  tho  mass  is  of  order  3 and  accelerations  aro 
also  of  order  p. 

Tho  toms  of  order  p yield  tho  following  conditions; 

§o  ® ° , 

2pg  J phdA  ■ M^pg  , 

J 

xphdA  - 0 , 


(1.15) 

(1.16) 
(1.17) 


(1.18) 


/ U|iflt-a0^lx)]t  dA  • 0 , 


(1.19) 


i’ 

A Cx^lt-8o^lx)^  cU  " 0 ' 


(1.20)  | Cy(^lt-s0/lx)!  clA  - o . 

The  symbol  ( occurring  horo  means  that  tho  Jump  in  the 
quantity  in  brackets  on  going  from  tho  positivo  tc  the  nogativo 
sido  of  tho  projoctod  aroa  A of  tho  ship’s  hull  is  to  bo  taken, 

Tho  variables  of  integration  are  x and  y.  The  oquation  (l*lf?) 
states  that  tho  term  of  order  zero  in  tho  speed  is  a constant,  and 
honco  tho  motion  in  tho  x-dirocticn  is  a small  oscillation  rela- 
tive to  a motion  with  unifemn  volccity.  Equation  (1.16)  is  an 
expression  of  tho  law  of  Archimodes*  tho  rest  position  of  equi- 
librium must  bo  such  that  tho  weight  of  tho  ship  just  equals  tho 
weight  of  the  water  it  displaces*  Equation  (1.17)  oxprossos 
another  law  of  equilibrium  of  a floating  body,  i.o,  that  the  cantor 
of  buoyancy  should  be  on  tho  same  vortical  lino  as  the  center  of 
gravity  of  tho  ship*  Tho  remaining  throo  equations  (1*18,  (1*19), 
and  (1*20)  in  the  group  serve  to  dotormino  tho  displacements 
©11#  ©21*  andOi,  which  occur  in  tho  boundary  condition  (1*12)  for 
tho  velocity  potential  fi i*  As  wo  havo  already  romarkod,  the  velo- 
city petontial  can  bo  written  in  tho  form  ( 1 * lip ) as  a linear 
combination  of  harmonic  functions  with  those  unknown  and  timo- 
dopondont  displacomonts  as  coefficients J insortion  in  oquations 


(ni8),  (1.19),  and  (1.20)  cloarly  leads  to  a coupled  system  of 

ordinary  differential  equations  with  constant  coefficients  for 

and  a^,  which  is  of  second  order  in  en#e21*  811(1  of 

first  ordor  in  fct,  (though  also  of  second  order  in  the  angular 

t • 

displacement  ^ ^ (*)£*»)•  7110  coefficients  of  those  difforon' 

tial  equations  are,  of  course,  obtained  in  terms  of  integrals  over 
A which  involve  tho  known  functions  ^^(x,y,0^.,t )•  If  the  spued 

s * const,  (which  occurs  in  (1.1k))  is  known,  it  follows  that 

o 

tho  systom  of  differential  equations  for  •id*)*  and  "L(t) 

would  yiold  those  displacements  uniquely  once  proper  initial 
conditions  aro  proscribed.  Wo  shall  see  in  a moment  that  sq  is 
fixed  by  a condition  that  is  independent  of  all  tho  unknown  dis- 
placements - in  fact,  it  depends  only  on  the  propoller  thrust  T 
and  tho  shape  of  tho  hull  - and  consequently  wo  have  obtained  a 
result  that  is  at  first  sight  rather  surprising!  tho  motion  of 
tho  wator,  which  is  fixed  solely  by  is  entirely  independent 
of  tho  pitching  displacement  ^^(t),  the  heave  yQ(t),  and  tho 
surge  s^(t),  i.o.  of  all  displacements  in  the  vortical  plane 
except  tho  constant  forward  spoed  s(_.  A libtlo  roflecticn, 
howovor,  makes  this  result  quite  plauslblo:  Our  tho  ry  is  basod 
on  tho  assumption  that  the  ship  is  a thin  disk  whero  thickness 
is  a quantity  cf  first  order  disposed  vertically  in  the  wator, 
hence  only  the  finite  displacements  of  tho  disk  parallel  to  this 
vortical  plane  - create  oscillations  in  tho  water  that  are  of 
socond  order  at  least.  On  tho  other  hand,  displacements  of  first 


order  of  tho  disk  at  right  anglos  to  itself  will  croato  motions  in 

tho  water  that  are  also  of  first  order#  Ono  might  sook  to  doscribo 

the  situation  crudsly  in  tho  following  fashion*  Imagine  a knifo 

blade  hold  vertically  in  tho  wator.  Up-and-down  motions  of  tho 

knifo  evidently  produce  motions  of  tho  water  which  aro  of  a quito 

different  ordor  of  magnitude  from  motions  produced  by  displacements 

of  tho  knifo  porpondlcular  to  its  blado.  Stress  is  laid  on  this 

phenomenon  horo  bocauso  it  holps  to  promoto  understanding  of  othor 

. * 

occurroncos  to  be  doscribod  lator* 

Tho  terms  of  second  ordor  in  (3  yield,  finally,  tho  following 
conditions* 


(1.21) 

(1.22) 


p j [hx(^Lt  * V*lx)]-  ^ + T» 


-2pg  i (y^  ♦ )hdx 


J (v-v! 


(1.23)  I31P31  * -2pg© 21  * (y-y^)hdA  - 2pgy1  J xhdx 
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- 2/g©3i  J x%idx  + IT  + p 1 IxlY"^-y^)^x3[/lt-so/lx)^dA 


We  no to  that  integrals  ever  tho  projootod  water-lino  L of  tho  ship 
in  its  equilibrium  position  occur  in  addition  to  intograls  over  tho 
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vortloal  projection  A of  tho  ontiro  hull*  The  quantity  X31  nr isos 

from  the  rolation  I *pi31  for  tho  moment  of  inertia  I of  the  ship 

with  rospoct  to  an  axis  through  its  c*g*  parallel  to  the  s' -axis* 

Tho  equation  (1*21)  determines  tho  surge  s^,  and  also  the  speed 

sQ  (or,  if  one  wishos,  the  thrust  T is  determined  if  sQ  is  assumed 

to  bo  given)*  Furthermore,  tho  speed  sc  is  fixed  solely  by  T and 

the  geomotry  of  the  ship's  hull*  This  can  bo  soon,  with  roforonco 

to  (1*14)  and  tho  discussion  that  aooompanies  it,  in  the  following 

way*  If  a>^,  -9'  , and  aro  constants,  then  they  must,  as  one 

could  show,  be  identically  soroj  hence  the  term  so^  in  ( 1* ll|. ) is 

the  only  term  in  that  is  independent  of  t*  It  therefore 

determines  T upon  insertion  of  ^ in  (1*21)*  This  term,  howovor, 

is  obtained  by  determining  the  harmonic  function  ^ as  a solution 

of  the  boundary  problem  for  in  tho  special  case  in  which 

* 

* ®21  * ^1  9 °*  h0nco,  as  ono  sees  from  (1,12)  and  (1*13), 

♦2  Is  fixed  by  sq  and  h^  alone*  In  fact,  the  relation  between  s^ 
and  T is  exactly  the  same  rolation  as  was  obtained  by  Micholl* 

(It  will  bo  written  down  lator* ) In  other  words,  tho  wavo  resis- 
tance is  now  soon  to  depohd  only  on  tho  basic  motion  with  uniform 
speed  of  the  ship,  and  not  at  all  on  its  small  oscillations  relative 
to  that  motioh*  If,  then,  offeots  on  tho  wavo  resistance  duo  to 
tho  oscillation  of  the  ship  are  to  bo  obtained  from  the  theory, 
it  will  bo  nooessary  to  take  acbount  of  higher  order  terms  in  order 
to  calculate  them*  Onoe  the  thrust  T has  boon  determined  tho 
equations  (1*22)  and  (1*23)  form  0 couplod  system  for  tho  determin- 
ation of  y^  and  ®ino©  ^ tuld  ®n  havo  presumably  boon 
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determined  previously*  However#  it  is  not  quite  correot  to  say  that 
the  surge  s^»  the  heave  y^#  and  the  pitching  osoillation  Qro 

not  oouplod  with  the  roll,  yaw  and  sway  sinoo  the  latter  quantities 
enter  into  the  determination  of  Thus  our  system  is  one  in 

which  there  is  a great  deal  of  cross-ooupling.  It  might  also  be 
noted  that  tho  trim,  i.e.  the  constant  values  of  y^  and  about 

which  tho  oscillations  occur  are  determined  from  (1*22)  and  (1.23) 
by  the  time -independent  terms  in  those  equations  — including  for 
example,  tho  moment  IT  of  the  thrust  about  the  c.g. 

We  havo  now  given  the  comploto  formulation  of  our  problem# 
except  for  initial  conditions  and  conditions  at  ® , Boforo  saying 
anything  about  motheds  for  finding  eonorotc  solutions  in  specific 
cases,  it  has  somo  point  to  mention  a number  of  conclusions,  in 
addition  to  those  already  givon,  which  can  bo  inforrod  from  our 
oquations  without  solving  them.  Consider,  for  cxamplo,  tho  equa- 
tions (1.22)  and  (1.23)  for  the  heavo  y^  and  tho  pitching  oscilla- 
tion ©-,#  and  make  tho  assumption  that  tho  integral  / xhdx  ■ 0 
J L 

(which  means  that  tho  horizontal  soction  of  tho  ship  at  tho  water 
lino  has  its  c.g.  on  the  same  vortical  as  that  of  tho  whole  ship). 

If  this  condition  is  satisfied  it  is  immodiatoly  soen  that  tho 
oscillations  and  y^  arc  not  coupled.  Furthermore#  those 

equations  are  seen  to  have  tho  form 
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(1.25) 

(1.26) 

with 

(1.27) 

(1.28) 

It  follows  immediately  that  rosonanco  is  pcssiblo  If  p(t)  has  a 
harmonic  component  of  the  form  A cos(X^t+B)  or  q(t)  a component  of 
the  form  A cosfXgt+B);  In  other  words,  ono  could  oxpoct  exceptionally 
heavy  oscillations • if  the  speed  of  tho  ship  and  the  sea  way  wore 
to  bo  such  as  to  load  to  forced  oscillations  having  frequencies 
close  to  those  valuos.  One  obsorvos  also  that  the so  resonant 
frequencies  can  bo  computod  without  reference  to  tho  motion  of  tho 
soa  or  tho  ship!  tho  quantitios  dopond  only  on  tho  shape  of 

the  hull.* 

In  spite  of  tho  fact  that  the  linear  thoory  presented  hero  must 
be  used  with  caution  in  relation  to  the  actual  practical  problems 

*Tho  equation  (1,27)  enn  bo  interpreted  in  tho  following  way;  it 
furnishes  tho  froquoncy  of  free  vibration  of  a system  with  one 
dogroo  of  freedom  in  which  tho  rostering  force  is  proportional  to 
tho  weight  of  water  displaced  by  a cylindor  of  cross*socti on  area 

2 J hdx  whon  it  is  immersed  vertically  in  water  to  a depth  y, • 

1 1 


h + xi?i  * p(t) 


®31  + X|  3i  * q(t) 


X?  = 


2pg  J hdx 

L 

M, 


2pg(  I (y-y^)hdA  + J x2hdx] 
A L 

hi 


V 
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concerning  ships  in  motion,  it  still  seoms  likely  that  such 

resonant  frequencies  would  be  significant  if  hqr  happened  to 

occur  in  the  terms  p(t)  or  q(t)  with  appreciable  amplitudes. 

Supposo,  for  instance  that  the  ship  is  moving  in  a sea-way  that 

consists  of  a single  train  of  simple  harmonic  progressing  plane 

waves  with  circular  frequency  <T  which  have  their  crests  at  right 

angles  to  tho  course  of  the  ship.  If  the  speed  of  the  ship  is  so 

one  finds  that  tho  circular  excitation  frequency  of  tho  disturbances 

caused  by  such  waves,  as  viewed  from  the  moving  coordinate  system 

s ,d2 

(x,y,z)  that  is  usod  in  tho  discussion  here,  is  (T  + —y—  , since 

O 

62 

is  2n  timos  tho  reciprocal  of  tho  wave  length  of  tho  wave  train# 

O 

Thus  if  X^  and  Xg  should  happen  to  lie  noar  this  value,  a heavy 

oscillation  night  be  expected.  One  can  also  soo  that  a change  of 

course  to  ono  quartering  tho  waves  at  angle  y would  load  to  a 

*2 

circular  excitation  frequency  f + bq  ers  y ' ^ and  naturally 

this  would  have  on  effect  on  tho  amplitudes  of  tho  rospenso# 

It  has  already  boon  stated  that  tho  work  presented  hero  is 
related  to  work  publishod  by  Haskind  [1]  in  1946*  and  it  was 
indicated  that  tho  two  theories  diffor  in  serno  rospocts.  Wo  havo 
not  mado  a comparison  of  the  two  theories  in  tho  gonoral  case, 
which  would  not  be  oasy  to  do,  b t it  is  possible  to  make  a compari- 
son rathor  easily  in  the  special  caso  treated  by  Haskind  in  dotail# 
This  is  the  special  case  troatod  in  tho  second  of  his  two  papers 
in  which  tho  ship  is  assumed  to  oscillato  only  in  the  vertical 
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plane  - as  would  be  possible  If  the  sea-way  oonslsted  in  trains  of 
plane  waves  all  having  their  crests  at  right  angles  to  the  course  of 
the  ship*  Thus  only  the  quantities  s^t),  y^ft),  and  ©^i^*  in 
our  notation,  would  figure  in  fixing  the  motion  of  the  ship* 

Haskind  treats  only  the  displacements  y^(t)  and  (which  are 

denoted  in  his  paper  by  £(t)  and  ty(t)),  for  which  he  finds  differ- 
ential equations  of  second  orderj  but  these  equations  are  not  the 
same  as  the  corresponding  equations  (1,22),  (1,23)  above.  One 

observes  that  (1,22)  contains  as  its  only  derivative  the  second 

• # 

derivative  y^  and  (1,23)  contains  as  its  sole  derivative  a term 
v 

with  other  words  there  are  no  first  derivative  terms 

at  all,  and  the  coupling  arises  solely  through  the  undifferentiated 
terms.  Haskind* s equations  are  quite  different  since  first  and 
second  derivatives  of  both  dependent  functions  occur  in  both  of 
the  two  equations j thus  Haskind,  on  the  basis  of  his  theory, 
can  speak,  for  example,  of  damping  terms,  while  the  theory  pre- 
sented here  does  not  yiold  any  such  terms.  The  authors  feci  that 
there  should  not  be  any  damping  terms  of  this  order  for  the  follow- 
ing reasons.  In  the  absence  of  frictional  resistances,  the  only 
way  in  which  onergy  can  be  dissipated  is  through  the  transport  of 
energy  to  infinity  by  means  of  out-going  progressing  waves.  How- 
ever, we  have  already  given  what  scorn  to  us  to  be  valid  reasons 
for  believing  that  the  oscillations  that  consist  solely  in  dis- 
placements parallol  to  the  vertical  piano  produoe  waves  in  the 
water  with  amplitudes  that  are  of  higher  order  than  thoso  considered 
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in  the  first  approximation.  Thus  no  such  dissipation  of  energy 

would  occur.  In  any  case,  our  theory  has  this  fact  as  one  of  its 

consequences.  Of  course,  it  does  not  matter  too  much  if  some 

terms  of  higher  order  are  included  in  a perturbation  theory,  at 

least  if  all  the  terms  of  a certain  given  order  are  really  present: 

at  most,  one  might  be  deceived  in  giving  too  much  significance  to 

the  higher  order  terms.  Haskind  also  says,  however,  and  we  quote 

from  the  translation  of  his  paper  (see  page  59)  * “'Thus,  for  a 

ship  symmetric  with  respect  to  its  midship  section  ........  only 

in  the  absence  of  translatory  motion,  i.e.  for  SQ  38  0,  are  the 

heaving  and  pitching  oscillations  independent."  This  statement 

does  not  hold  in  our  versior^of  the  theory.  As  one  sees  from 

(1.22)  and  (1.23)  coupling  occurs  if,  and  only  if  I xhdx  + 0, 

L 

whether  S vanishes  or  not.  In  addition,  Kaskind  obtains  no 
o 

resonant  frequencies  in  these  displacements  because  of  the  presence 
of  first-derivative  terms  in  his  equations;  the  authors  feel  that 
such  resonant  frequencies  may  well  be  an  important  feature  of  the 
problem. 

We  turn  next  to  a brief  discussion  of  methods  of  solving  the 
problems  formulated  above.  The  difficulties  are  for  the  most  part 
concentrated  in  the  problem  of  determining  the  first  approximation 
/^(x,y,z,t)  to  the  velocity  potential.  The  discussion  above  assumed 
that  ^ had  in  some  way  been  determined  in  the  form  (1.14)  by 
solving  the  boundary  value  problem  posed  by  (1.11),  (1.12),  (1.13), 
and  appropriate  conditions  at  the  time  ts0  and  at  ® • In  general, 
an  explicit  solution  of  the  problem  for  - in  terms  of  an  integral 
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representation,  say  - seems  out  of  the  question*  In  fact,  as  soon 
as  rolling  or  yawing  motions  occur,  explicit  solutions  are  unlikely 
to  be  found*  The  best  that  the  authors  have  been  able  to  do  so 
far  in  such  cases  has  been  to  formulate  an  integral  equation  for 
the  values  of  ^ over  the  vertical  projection  A of  the  ship's  hull; 
this  method  of  attack,  which  looks  possible  and  somewhat  hopeful 

I 

for  numerical  purposes  since  the  motion  of  the  ship  requires  the 
knowledge  of  ^ only  ovor  the  area  A,  is  under  investigation* 
However,  If  the  motion  of  the  ship  is  confined  to  a vertical  plane, 
so  that  6J1  = ©^  = ©*21  = 0,  it  is  possible  to  solve  the  problems 

explicitly.  This  can  be  seen  with  reference  to  the  boundary 
conditions  (1.12)  and  (1.13)  which  in  this  case  are  identical  with 
those  of  the  classical  theory  of  Michell  and  Hsv clock,  and  hence 
permit  an  explicit  solution  for  ft,  which  is  given  later  on  in 
section  I4..  Aftor  ^ is  determined,  it  can  be  inserted  in  (1.21), 
(1.22),  and  (1.23)  to  find  the  forward  speed  S0  corresponding  to 
the  thrust  T,  the  two  quan titles  fixing  the  trim,  and  the  surge, 
pitching,  and  heaving  oscillations.  In  all,  six  quantities 
fixing  the  motion  of  the  ship  are  determined. 

The  theory  developed  in  this  report  is  very  general,  and  it' 
therefore  could  bo  applied  to  the  study  of  a wide  variety  of 
different  problems*  For  example,  tho  stability  of  the  oscillations 
of  a ship  could  be  investigated  on  a rational-dynamical  basis, 
rather  than  as  at  present  by  assuming  the  water  to  remain  at  rest 
where  the  ship  oscillates*  It  would  be  possible  in  principle 
to  inveatigate  theoretically  how  a ship  would  move  with  a given 
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rudder  setting,  and  find  the  turning  radius,  angle  of  heel,  etc. 
The  problem  of  stabilisation  of  a ship  by  gyroscopes  or  other 
devices  could  be  attacked  in  a very  general  way:  the  dynamical 
equations  for  the  stabilizers  would  simply  be  included  in  the 
formulation  of  the  problem  together  with  the  forces  arising  from 
the  interactions  with  the  hull  of  the  3hip. 

In  Sec.  2 the  general  formulation  of  the’  problem  is  given; 
in  Sec.  3 the  details  of  the  linearization  process  are  carried 
out;  and  in  Sec.  1+  a solution  of  the  problem  is  given  for  the  case 
of  purely  vertical  motion,  including  a verification  of  the  fact 
that  the  wave  resistance  is  given  by  the  same  formula  as  was 
found  by  Michell. 

2.  General  formulation  of  the  problem. 

Wo  derive  here  the  basic  theory  for  the  motion  of  a floating 
rigid  body  through  water  of  Infinite  depth.  The  water  is  assumed 
to  be  in  motion  as  the  result  of  the  motion  of  the  rigid  body, 
and  also  because  of  disturbances  at  00 ; the  combined  system  con- 
sisting of  the  rigid  body  and  the  water  is  to  be  treated  as  an 
interaction  in  which  the  motion  of  the  rigid  body,  for  example, 
is  determined  through  the  pressure  forces  exerted  by  the  water  on 
its  surface.  We  assume  that  a velocity  potential  exists.  Since 
we  deal  with  a moving  rigid  body  it  is  convenient  to  refer  the 
motion  to  various  types  of  moving  coordinate  systems  as  well  as 
to  a fixed  coordinate  system.  The  fixed  coordinate  system  is 
denoted  by  0-X,Y,Z.  The  X,Z-plano  is  in  the  equilibrium  position 
of  the  free  surface  of  the  water,  and  the  Y-axis  is  positive 
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upwards.  The  first  of  the  two  moving  coordinate  systems  we  use 

(the  second  will  be  introduced  later)  is  denoted  by  o-x,y,z,  and 

is  specified  by  Fig.  1.1.  The  x,z-plane  coincides  with  the  X,  Z- 

plano  (i.e.  it  lies  in  the  undisturbed  free  surface),  the  y-axis 

is  vertically  upward  and  contains  the  center  of  gravity  of  the 

ship.  The  x-axis  has  always  the  direction  of  the  horizontal 

component  of  the  velocity  of  the  center  of  gravity  of  the  ship. 

(If  we  define  the  course  of  the  ship  as  the  vertical  projection 

of  the  path  of  its  center  of  gravity  on  the  X, Z^plane,  then  our 

convention  about  the  x-axis  means  that  this  axis  is  taken  tangent 

to  the  ship’s  course.)  Thus  if  R * (X  ,Y  ,Z  ) is  the  position 

c c c c 

vector  of  the  center  of  gravity  of  the  ship  relative  to  tho  fixed 

*•  ... 

coordinate  system  and  hence  R 3 (X  ,Y  ,Z  ) is  the  velocity  of 

c c c c 

the  c.g.,  it  follows  that  the  x-axis  has  the  direction  of  the 
vector  ft  given  by 


(2.1) 


u = X I + Z J 
c c 


with  I and  J unit  vectors  along  the  X-axis  and  the  2-axis.  If 
l is  a unit  vector  along  the  x-axis  we  may  write 


(2.2) 


s(t)i  = u 
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thus  introducing  the  speed  s(t)  of  the  ship.  For  later  purposes 
we  also  introduce  the  angular  velocity  vector  u of  the  moving 
coordinate  system: 


t 

* 


and  the  angle  a (cf,  Pig,  1,1)  by 


t 

(2,i|.)  a(t ) * J w(T}dT  , 

0 

Tlie  equations  cf  transformation  from  one  coordinate  system 
to  the  other  are 


(X  * x cos  a+z  sin  a+Xc  i x 

X = y jy 

Z **  -x  sin  a+t  cos  a+Z  j z 

c 


* (X-X  ) cos  a-( Z-Z  ) sin  a 
c c 

= Y 

= (X-X  ) sin  a+(Z-Z  ) ccs  a < 
c c 


L H 


t 1 
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By  f(x,  Y,Z;t)  we  denote  the  velocity  potential  and  write 

(2,6)  $(X,Y,Zjt)  * f(x  ccs  a+z  sin  a+XQ,y,-x  sin  a+z  cos  a+ZQ;t) 

5 4 ( x,  y , z j t ) . 

In  addition  to  the  transformation  formulas  for  the  coordinates,  we 
also  need  the  formulas  for  the  transformation  cf  various  derivatives. 
One  finds  without  difficulty  the  following  formulas: 


(2.7) 


4X  cos  a+4z  cos  a 


<$Y 


Ty 

- A 8in  a+4„  cos  a 
X z 


4* . ; 


' , X 


v 


It  is  clear  that  grnd2  $(X,Y,Z;t)  * grad2  4(x,y,z;t)  and  that  4 
is  a harmonic  function  in  x,y, z since  § is  harmonic  in  X,Y,Z.  To 
calculate  |i3  a little  more  difficult;  the  result  is 

(2,6)  §£*  -(s  +t*iz)4x  +yxl2  + . 

(To  verify  this  formula,  one  uses  a <i*xxt+^yyt+^zzt+^t  ‘'la0 

relations  (2*5)  together  with  s cos  a * X_,  s sin  a * -Z  * ) 

c c 

The  last  two  sets  of  equations  make  it  possible  to  express  Ber- 
noulli ' s law  in  terms  of  <|i(x,y,zt); 

(2.9)  ^ + gy  + ^(grad  4)2  + (s+uz)4x  -wx<J>z  - <|>t  * 0 . 

Suppose  now  that  P(X,Y,Z;t)  = ‘0,  is  a boundary  surface 
(fixed  or  moving)  and  set 

(2.10)  F(x  cos  a+...,y,-x  sin  a+...;t)  2 f(x,y,z;t)  , 

so  that  f(x,y,z;t)  s o is  the  equation  of  the  boundary  surface 
relative  to  the  moving  coordinate  system.  The  kinematic  con- 
dition to  be  satisfied  on  such  a boundary  surface  is  that  the 

dF 

particle  derivative  vanishes,  and  this  leads  to  the  boundary 
condition 

(2.11)  4xfx  + + 4zrz  - *(s+wz)fx  + wxfz  + ft 


relative  to  the  moving  coordinate  system  upon  using  the  appropriate 
transformation  formulas.  In  particular,  if  y - \(x,  Z|t)  * 0 is 
the  equation  of  the  free  surface  of  the  water,  the  appropriate 
. kinematic  condition  is 

(2.12)  -4xnx  + ^y  “ ^znz  “ (s+<*>2)nx  - o*\.z 

to  be  satisfied  for  y = ^ • (The  dynamic  free  surface  condition 
is  of  course  obtained  for  y s H from  (2.9)  by  setting  p = 0.) 

We  turn  next  to  the  derivation  of  the  appropriate  conditions, 

both  kinematic  and  dynamic,  on  the  ship*s  hull.  To  this  end  it 

is  convenient  to  introduce  another  moving  coordinate  system 

o*  - x^y’jZ*  which  is  rigidly  attached  to  the  ship.  It  is 

assumed  that  the  hull  of  the  ship  has  a vertical  plane  of  symmetry 

(which  also  contains  the  center  of  gravity  of  the  ship);  we  locate 

the  x'jy'-plane  in  it  (cf.  Pig  1.2)  and  suppose  that  the  y»-axis 

contains  the  center  of  gravity.  The  o’-x^y*,2'  system,  like 

the  other  moving  system,  is  supposed  to  coincide  with  the  fixed 

system  in  the  rest  position  of  equilibrium  at  t=0.  The  center  of 

gravity  of  the  ship  will  thus  be  located  at  a definite  point  on 

the  y’-axis,  say  at  distance  y£  from  the  origin  0» : in  other 

words,  the  system  of  coordinates  attached  rigidly  to  the  ship  is 

such  that  the  center  of  gravity  has  the  coordinate  (0,y>,0). 

c 

In  the  present  section  we  do  not  wish  in  general  to  carry 
out  linearizations.  However,  since  we  shall  in  the  end  deal  only 
with  motions  which  involve  small  oscillations  of  the  ship  relative 
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to  the  first  moving  coordinate  system  o-x,y,z,  it  is  convenient 
and  saves  time  and  space  to  suppose  even  at  this  point  that  the 

angular  displacement  of  the  ship  relative  to  the  o-x,y,z  system 

*> 

is  so  small  that  it  con  be  treated  as  a vector  If: 

(2.13)  0 * ©!?  + ©2J  + ©3k  . 

The  transformation  formulas,  correct  up  to'  the  first  order  terms 

*» 

in  the  components  0^  of  0,  are  then  given  by: 

X*  » X + 03(y-yc)  - »22 

(2.14)  y*  = y * (yc-yj)  + ©jz  - #3x 
*«  = z + ©2x  - ®i<yyc) 

with  yQ  of  course  representing  the  y-coordinnte  of  the  center  of 
gravity  in  the  unprimed  system*  It  is  assumed  tyat  y - y»  is  a 
small  quantity  of  the  same  order  as  the  quantities  ©i  and  only 
linear  terms  in  this  quantity  have  been  retained*  (The  verifi- 
cation of  (2.1 4)  is  easily  carried  out  by  making  use  of  the  vector- 
*•  #•  #►  *■ 
product  formula  0 “ 0 X r,  for  the  small  displacement  6 of  a rigid 

body  under  a small  rotation  0. ) 

The  equation  of  the  hull  of  the  ship  (assumed  to  be  symmetrical 
with  respect  to  the  x* ,y* •plane ) is  now  supposed  given  relative  to 
the  primed  system  of  coordinates  in  the  form: 


(2.15) 


> 
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z*  * - £(x' ,y’ ) , z»  ^ 0 . 

The  equation  of  the  hull  relative  to  the  0-x,y, z-systom  can  be 
written  in  the  form 

(2.16)  z+e2x-  ^(y-y’)  - S(x,y)  + [©2z-©3(y-y* )]?x(x>y) 


+X(yc-yJ)-V+e3x]S(x,y)  = 0 » z*  > 0 * 


when  higher  order  terms  in  (y  -y*)  and  0,  arc  neglected*  The  left 
hand  side  of  this  equation  could  now  be  inserted  for  f in  (2.11) 
to  yield  the  kinematic  boundary  condition  on  the  hull  of  the  ship, 
but  we  postpone  this  step  until  the  next  section. 

The  dynamical  conditions  on  the  ship’s  hull  are  obtained 

from  the  assumption  that  the  ship  is  a rigid  body  in  motion  under 

> *• 

the  action  of  the  propeller  thrust  T,  it3  weight  Mg j,  and  the 
pressure  p of  the  water  on  its  hull.  The  principle  of  the  motion 
of  the  center  of  gravity  yields  the  condition 

. * ♦»  f ► <►  #» 

(2.17)  M 3pai  + ycJ)  * J p n dS  + T • MgJ  . 

S 


By  n we  mean  the  inward  unit  normal  cn  the  hull, 
coordinate  system  o-x,y, z is  such  that  ^ * -«k 
so  that  (2.17)  can  be  written  in  tho  form 


Our  moving 

and  ^|*0, 
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(2.18)  Mai  - Ms  w k + MycJ  ■ J p n dS  + T -MgJ  , 

S 

with  p defined  by  (2.9).  The  law  of  conservation  of  angular 
momentum  is  taken  in  the  form: 

. ( ► ► *■  ► r ► ♦ •+•  -►  •*» 

(2.19)  j (R-Rc)  X(R-R0)dm  = J p(R-Rc)Xn  dS  + (yRc)XT  . 

M S 

The  crosses  all  indicate  vector  products.  By  R is  meant  the 
position  vector  of  the  element  of  mass  dm  relative  to  the  fixed 
coordinate  system.  Rc  fixes  tho  position  of  the  C.g.  and  R^, 
locates  the  point  of  application  of  the  propeller  thrust  T,  also 
relative  to  tho  fixed  coordinate  system.  Wo  introduce  r = (x, y,z) 
as  the  position  vector  of  any  point  in  the  ship  relative  to  the 
moving  coordinate  system  and  set 

(2.20)  q * r - ycj  , 

*• 

so  that  q is  a vector  from  the  c.g.  to  any  point  in  the  ship.  The 
relation 

• • • 

> > > > > 

(2.21)  R = Rc  + (u»  + 9)  x q 

holds,  since  + © is  the  angular  velocity  of  the  ship;  thus 

(2.21)  is  simply  the  statement  of  a basic  kinematic  property  of 

* 

rigid  bodies.  By  using  the  last  two  relations  the  dynamical 


condition  (2.19)  can  be  exprossod  in  terns  of  quantities  measured 
with  respect  to  the  moving  coordinate  system  o-x,y, z,  as  follows: 

(2.22)  |t  / (r*ycJ)  X [(«•>  + ©)  X (r-ycJ)]  dm 

* J P(r-yeJ)*n  dS  + (Rt  • Re>  H I . 
s 

We  have  now  derived  the  basic  equations  for  the  motion  of 
the  ship.  What  would  be  wanted  in  general  would  be  a velocity 
potential  <J>(x,y, z;t)  satisfying  (2.31)  on  the  hull  of  the  ship, 
conditions  (2.9)  (with  p * 0)  and  (2.12)  on  the  free  surface  of  the 
water;  and  conditions  (2.1?)  and  (2.22),  tdiich  involve  4 under 
integral  signs  through  the  pressure  p as  given  by  (2.9).  In 
addition,  there  would  be  initial  conditions  and  conditions  at  od 
to  be  .satisfied.  Detailed  consideration  of  these  conditions,  and 
the  complete  formulation  of  the  problem  of  determining  <j:(x,y,z;t ) 
under  various  conditions  will  be  postponed,  however  until  later  on 
since  we  wish  to  carry  out  a linearization  of  all  of  the  conditions 
formulated  here. 

3.  Linearization  by  a formal  perturbation  procedure. 

Because  of  the  complicated  nature  of  our  conditions,  it  seems 
wise  to  carry  cut  tho  linearization  by  a formal  development  in 
order  to  make  sure  that  all  toms  of  a given  order  are  retained; 
this  is  all  the  more  necessary  since  terns  of  different  orders 
must  be  considered.  The  linearization  carried  out  here  is  based 
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on  the  assumption  that  the  motion  of  the  water  relative  tc  the 
fixed  coordinate  system  is  a small  oscillation  about  the  rest 
position  of  equilibrium.  It  follows,  in  particular,  that  the 
elevation  of  the  free  surface  of  the  water  should  be  assumed  to  be 
small.  We  do  not,  however,  wish  to  consider  the  speed  of  the 
ship  with  respect  tc  the  fixed  coordinate  system  to  be  a small 
quantity:  it  should  rather  be  considered  a finite  quantity.  This 
brings  with  it  the  necessity  to  restrict  the  form  of  the  ship  so 
that  its  motion  through  the  water  does  net  cause  disturbances  so 
large  as  to  violate  our  basic  assumption;  in  other  words,  we  must 
assume  the  ship  to  have  the  form  of  a tbin  disk.  In  addition,  it 
is  clear  that  the  veloci  ty  of  such  a disk-like  ship  must  of  neces- 
sity4maintain  a direction  that  does  net  depart  too  much  from  the 
plane  of  the  thin  disk  if  only  small  disturbances  in  the  water  are 
to  be  created  as  a result  of  its  motion  with  finite  speed.  Tims 
we  assume  that  the  equation  of  the  ship’s  hull  is  given  by 

(3.1)  z'  * ph(x',y»)  , z’  > 0 , 

with  p a small  dimensionless  parameter,  so  that  the  ship  is  a 
thin  disk  symmetrical  with  respect  to  tho  x’ ,y’ -plane,  and  ph 
takes  the  place  of  £ in  (2.15)*  (It  might  be  noted  in  passing 
that  this  is  not  the  most  general  way  to  describe  the  shape  of  a 
disk  that  would  be  suitable  for  a linearization  of  the  type  carried 
out  here.)  We  have  already  assumed  in  the  preceding  section  that 
the  motion  of  the  ship  is  a small  oscillation  relative  to  the  moving 


coordinate  system  o-x,y,z  - an  assumption  that,  in  fact,  is  made 
necessary  by  our  basic  assumptions  concerning  the  linearization. 
It  seems  reasonable,  therefore,  tc  develop  all  of  our  basic 
quantities  (taken  as  functions  of  x,y,z;t)  in  powers  of  p,  as 


follows: 

(3*2) 

4(x,y,z;t;p) 

= pij  + p2^ 

(3*3) 

Yl(x,z;t;p) 

= 

(3.U 

s(t;p) 

= s + ps^  + 

13*5) 

U)  (t  jp ) 

= coo  + Po»i  + 

(3.6) 

^(tfP) 

= pfin  + p2e. 

(3.7) 

yc  * yc 

= <*!  ♦ P2y2 

12 


The  first  and  second  conditions  state  that  tho  velocity  potential 
and  the  surface  wave  amplitudes,  as  seen  from  the  moving  system, 
are  small  of  order  p«  The  speed  of  the  ship,  on  the  ether  hand, 
and  the  angular  velocity  of  its  c.g.  about  the  vertical  axis  of 
the  fixed  coordinate  system,  are  assumed  tc  be  of  cade  r zero,  (It 
will  turn  out,  however,  that  coq  * 0 - a not  unexpected  result.) 
The  relations  (3.6)  and  (3*7)  serve  tc  make  precise  our  previous 
assumption  that  the  motion  of  the  ship  is  a small  oscillation 
relative  to  the  system  o*x,y,z. 

We  must  now  insert  these  developments  in  the  conditions 
derived  in  the  previous  section.  The  free  surface  conditions  are 
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treated  first.  As  a preliminary  step  we  observe  that 

(3.8)  4x(x,n,z; t;p)  * PUlx(x,o,z;t)  + n.4lxy(x,o,z|t)  + ...] 

* P2^2x+^2xy+ 

+ 

* P*lx<x'  °*  )+P2tni*i^(x,  0,  z; t W2x(x,  0, z; t ) j 


with  similar  formulas  fcr  other  quantities  when  they  are  evaluated 
on  the  free  surface  y = X • Here  we  have  used  the  fact  that  X 
is  small  of  order  p and  have  developed  each  term  in  a Taylor  series. 
Consequently,  the  dynamic  free  surface  condition  fcr  y “ Y\  arising 
from  (2.9)  with  p 3 0 can  be  expressed  in  the  form 

g(?li1+P2»l2+  + ^(32[(grad  + •••] 

+ [s0+Ps1+...+z(qo+P*^  +.../] C P^ix+P2(^i^ixy+^2x)+* * * 1 

(3.9) 

- x(oQ+p<^+. . . ) [ P$lz+P2 (nL+i^ag)*. • • 1 

-[P^t+P^ity+^t  )♦...)  * 0 

and  this  condition  is  to  be  satisfied  for  y * 0*  In  fact,  as 
always  in  problems  of  small  oscillations  of  continuous  media,  the 
boundary  conditions  are  satisfied  in  general  at  the  equilibrium 
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position  of  the  boundaries.  Upon  equating  the  coefficient  of  the 
lowest  order  tern  to  zero  we  obtain  the  dynamical  free  surface 
condition 

(3*10)  S»t1+(so+ucz^lx  “ °cx<^lz"^lt  = 0 fcr  7 * 0 » 

and  it  is  clear  that  conditions  on  the  higher  order  terms  could 
also  be  obtained  if  desired.  In  a similar  fashion  the  kinematic 
free  surface  condition  can  be  derived  from  (2,12) j the  lowest  order 
term  in  {3  yields  this  condition  in  the  form: 

(3.11)  <*ly  - (s0+  Mcz)\lx  +u0x\lz+\t  * 0 for  y * 0 . 

We  turn  next  to  tho  derivation  of  the  linearized  boundary 
conditions  on  the  ship's  hull.  In  view  of  (3*6)  and  (3.7),  the 
transformation  formulas  (2.14)  can  be  put  in  the  form 

/**  * x+p«31(y-y*)  -P»21z 

(3.12)  jy*  * y-^y1+po11z-pe31x 

* z+P®21x“P®ll*y“ycJ 

when  terms  involving  second  and  higher  powers  of  (3  are  rejected. 

Consequently,  the  equation  (2.16)  of  the  ship's  hull,  up  to  terms 
2 

in  p , can  be  written  as  follows 
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Z+P®21X"P®11 ( y"yc ) -Phf x+(3©3i  (y-y£ ) -P®21z,  y-py1+p«11E-p®31x]  = 0, 


and,  upon  expanding  the  function  h,  the  equation  becomes 


(3.13)  z+P©21x-pei;L(y-yi  )-ph(x,y)+..«  * 0 , 


tho  dots  representing  higher  order  terns  in  p.  We  can  new  obtain 
the  kinematic  boundary  condition  fc  r the  hull  by  inserting  the 
left  hand  side  of  (3.I3)  for  the  function  f in  (2.11);  the  result 
is 


(3.34) 


Ulz  * “sc(®21-hx)+x  l+®21x’®ll(y’yc) 


when  the  terms  of  zero  and  first  order  only  are  taken  into  account. 
It  is  clear  that  these  conditions  are  to  be  satisfied  over  the 
domain  A of  the  x,y-plane  that  is  covered  by  the  projection  of  the 
hull  on  the  plane  when  the  3hip  is  in  the  rest  position  cf  equili- 
brium. As  was  mentioned  earlier,  it  turns  out  that  «=  0,  l.e. 
that  the  angilar  velocity  about  the  z-axis  of  the  c.g.  of  the  ship 
in  its  course  must  be  small  of  first  order,  or,  as  it  could  also 
be  put,  the  curvature  of  the  ship’s  course  must  be  small  since  the 
speed  in  the  course  is  finite.  The  quantity  s^t)  in  ( 3 *4 ) thus 
yields  the  oscillation  of  the  ship  relative  to  the  x-axis. 


It  should  also  be  noted  that  if  we  use  z - -ph(x,y)  we 
find,  corresponding  to  (3.H4.),  that 

Kz  - -ao(*2i+V  + <"i+e2i)x  - “n'y-yJ1  • 

This  means  that  A must  be  regarded  as  two  sided,  and  that  the 

last  equation  is  to  be  satisfied  on  the  side  of  A which  faces 

the  negative  z-axis.  The  last  equation  and  (3 *14)  implies  that 

p may  be  discontinuous  at  the  disk  A, 

The  next  step  in  the  procedure  is  to  substitute  the  develop 

ments  with  respect  to  p,  (3*2)  - (3»7)»  in  the  conditions  for  the 

ship*s  hull  given  by  (2.18)  and  (2.22).  Let  us  begin  with  the 

integral  / p n ds  which  appears  in  (2.18).  In  this  integral  3 
S 

is  the  immersed  surface  of  the  hull,  n is  the  inward  unit  normal 
to  this  surface  and  p is  the  pressure  on  it  which  is  to  be  calcu- 
lated from  (2.9).  With  respect  to  the  o-x,y,z  coordinate  system 
the  equations  of  the  symmetrical  halves  of  the  hull  are 

S1s-  z * H^x^tip)  = fx  + f2 

(3.15) 

V-  z * H2(x,y,t;p)  = -f1  + f2 


where 


(3.16) 


t ! = ph  + p C«31(y-y^V(e3ix+yi)hy3  + 0((3  } 

f2  * -P©21X  + pe11(y-y»)  + 0(p2)  . 


We  oan  now  write 


J p £ ds  * / p Vn  + / p n2di 


in  which  and  n2  are  given  by 


>•  ♦> 


n. 


H..  1+H.  j-k 
_ lx  lyJ 


,1+!'L+Hly 


n. 


Vl+H?.+H2 


2x  2y 


We  can  also  write 


J p n ds  * -pg  J y n ds  + J p^n  ds  = 

s s s 


* -pe  J 7 n ds  + y P1n1ds1  + / Px  n 


ds 


where  p^,  from  (2.9),  is 


(3.17)  **  -p[^(grad  /02  + (s+oz)^x  - xl*(z  - ^3* 
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If  is  the  hull  surface  below  the  xz-plane,  the  surface 
o 

area  SQ-S  is  of  order  p and  in  this  area  each  of  the  quantities 
y,  Hp  H2  is  of  order  0.  Hence 

- *■  ' + *»  +■  -1*2 
- J y n ds  » - J y n ds  + (i.+j)  0 (p-5)  + k 0(  p*) 

S S 

o 

Prom  the  divergence  theorem 

- J y n ds  * 

So 

where  V is  the  volume  bounded  by  SQ  and  the  xz-plane.  With  an 
accuracy  of  order  p^,  V is  given  by 

V ■ 2p  / hdA  - f p(y1+«3lx)dB  « 2p  J hdA  - 2p2  /(y^+xd^ )hdx 
A B A L 

Here  A is  the  projection  of  the  hull  on  the  vertical  plane  when  the 
hull  is  in  the  equilibrium  position,  B is  the  equilibrium  water  line 
area,  and  L is  the  projection  of  the  equilibrium  water  line  on  the 
x-axis. 

If  W^,W2  are  the  respective  projections  of  the  immersed 
surfaces  S^,  S2  on  the  xy-plane  we  have 


iVw 


<u 
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] Px  n d - i / j P1(x#y,H1it)HlxdW1-|  P^x^Hgjt )H2xdW2  j 

s V.W  w2  J 

+ jj  / P1(x#y*H1|t)HlydW1-  Pj^tx^jHgjtjHgydwJ 

lwi  W2 

" k|  / P1(x»y*H1it)dWi  “ / PxCx^yjH^tJdWg  J 

lw,  W-  ' 


Neither  nor  W2  is  equal  to  A*  Each  of  the  differences  W^-A, 
W2**A  is,  however,  an  area  of  order  (3«  Prom  this  and  the  fact  that 
each  of  p,  H^,  H^y,  H2x,  H2y  is  of  order  P;  it  follows  that 

[Pi  ( x,  y,  H^*  t ( x,  y,  K2;  t )H2x]  dA+0  ( p3 ) j 

+ j|/[P1(x,y,H1;t)Hly-p1(x,y,H2;t)H2y]dA+0(p3)| 

- ^|j[p1(x,y,H1jt)-p1(x,y,H2;t)]dA+0(p2)  } 


f Pi  n"ds  * i{/ 
s i a 


It  was  pointed  out  above  that  / may  be  dlscondinuous  on  A» 
Hence  from  (3.17),  (3.2),  (3.4) 


(3.19) 


PiU^H^t)  * pP^it,s0^ixJ+  + 0 
P!(*#y,H2|t)  * PP^it"so^lx*"  + °^2)  • 


Here  the  + and  • superscripts  denote  values  at  the  positive  and 
negative  sides  of  the  disk  A whose  positive  side  is  regarded  as 
the  side  which  faces  the  positive  z-axis.  If  we  substitute  the 
developments  of  Hlx,  Hly,  H2x,  H2y,  and  (3.19)  in  (3.18),  then 
collect  the  previous  results,  we  find 

fa  n ds  * i{pP2/t(Ve2i,('*n*ao'fix)++(Ve2i)('fit-3o^ajtridA+0(P3)i 

3 t A 


2pgpJhdA  - 2pgp2J ( y1+x©^1 )hdx 


(3.20)  +PP2/t  ^ V®H ) (^lt“so^lx)  +(hy“eil)(^lt*so,6lx)  3dA*0(P 

^ A 

“lc  |ppj[  (^lt“So^lx^  * ^It-Vlx*  )dA+°(P^)^  . 


The  integral  /p(^-y^)*  £ ds  which  appears  in  (2.22)  can  be 


Jp(r-ycJ)x  & ds  * -pg  J y(r-yc"?)x  n ds 


+ J P1(r-ycJ)X  nx  dsx 


r -4  *4.  -4 

+ / P. (r-y_ j)x  n~  dS„ 


If  we  use  the  same  procedure  as  wa3  used  above  for  the  expansion 

of  / p n ds  we  find 
S 


► *- 


f p(r-ycJ)x  n ds 

s 


* -i|pp/[(y-yc)(^it-so^ix)+’(y"yc)(^it“so^ixr3dA+0(p2)} 

+j^p/[x(^lt-so^lx)+-x(/lt-3o^lx)"]dA+0(p2)j 

2pgP.^chdA-2ogp2  «3i/(y-y')h  dA 

-ZpgP^Jxhdx-apgP2#  ^/x2hdx 

L L 

+K^,pP‘"JT[x(hv+©1i)  (^it“a0^lx3  +x^hv”®H3  ^lt”so^lx3  3dil 

ii  * v 

-pP2/t(7-ycHV821)('(lt-so'flx'++(7'7c)(V821)(,fn*3e^lx)‘>J'‘ 

ifc 

+ 0(P3) 

•*'“  ' 2 

We  now  assume  that  the  propeller  thrust  T is  cf  order  p 
and  is  directed  parallel  to  the  x'oxi^t  that  is 


I S P2  I il 
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where  1*  is  tho  unit  vector  along  the  x'-axis*  We  also  assume  that 

#- 

T is  applied  at  a point  In  the  longitudinal  plane  of  symmetry  of 
the  ship  Jl  units  below  the  center  of  mass*  It  then  follows  that 


(3.22) 


T * p2T3  • + 0(p3) 


and 


(3.23) 


» ► *• 

(R  - R )x  T = • JL  JxT 


fl  «2n 


The  mass  of  the  ship  is  of  order  p*  If  wo  write  M = M, P and 
expand  tho  left  hand  side  of  (2.18)  in  powers  of  p it  becomes 


(3.24) 


. [i;1pso+I41p2s1+0(  p3 ) ] + J (M1p2'y1+o  ( p3 ) ) -k[  0(  p3 ) ] 

t *■  *■  *■ 

»Jpnds+T  - IL  Pg J • 


The  expansion  of  tho  left  hand  side  of  (2.22)  gives 


i 1 0(  p2 ) ] + j [ o ( p2 ) ] +Tc[  i31P2e'31+o  ( p3 ) 3 
P(p-ycJ)*  n ds  + (RT-Rc)jr  T 


(3.25) 


4.8 


where  is  the  moment  of  inertia  cf.  the  ship  about  the  axis 

which  is  perpendicular  to  the  longitudinal  plane  cf  symmetry  cf 
the  ship  and  which  passes  through  the  center  cf  mass. 

If  we  replace  tho  pressure  integrals  and  thrust  terms  in  the 
last  two  equations  by  (3*20),  (3.21),  (3*22),  (3.23)»  and  then 
equate  the  coefficients  cf  like  powers  of  {5  in  (3.24)  and  (3*25) 
we  obtain  the  following  linearized  equations  of  motion  of  the  ship. 
From  the  first  order  terms  we  find 

(3.26)  s0  * 0 

(3.27)  2pg  / |3hdA  * M-pg 

A A 

(3.28)  . /xphdA  = 0 

A 

(3.29)  /t^it’so^lx,+"(^lt“so^lx,’]dA  * 0 

(3.30)  J fx(^lt"ao^lx)+-x^lt'so^lxridA  = 0 

A , 

(3.31)  / [(y-y')(/lt-s0/1?r)+-(y.y»)(/lt.sc/lx)“]dA  * 0 

or  by  (3.29) 

(3.32)  / ty (^lt-s0^ix ) *-7 <*it 'so*lx * " J dA  “ 0 * 

A 


''"tT' 


V 


Jk* 


Prom  the  second  order  terms  we  find 

Vi  - p|t(V#2i)(^t-Bo'fi*)++(V®ai,<'fit'so|l'i*)‘)aA+T 

(3.33) 

= pjlh x(<*lt'so^lx)++hx(^lfso^lx>’14A+T 

Vi  * -2re/(yi^e3i)hdx 

Lt 

+p/[(hy+eu){/lt-s0^lx)++(hy-e11)(/lt-s0^lx)  ]dA 

(3.34)  3 -2pg/(y1+x93l)hdx 

L 

+p/[hy(^lt’so^lx)++hy(^lt‘so^lxr3dA 
n J 

r3i®3i  “ -2Ps®3i/<y*yJ)MA-2Pw1/*hd* 

-2pg  0^fx2hdx+  J T 

*p/(x<V^  (^lt*,<A*)++x(  V*«  ’ ^lf*<Ax>‘laA 

-pit  ly-yj ) < V*2i ) ('*xt*,o'fix)++(3r"yc 1 ( V®21  ‘ A 

A 

or  by  (3.30),  (3*31) 


1 
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(3.35) 


I31®31 


-2pg«31/(y-y^)hdA  - 2pgy1/xhdx 

+ Jt  T 

+p/[*hy-(y-y'»g  t (/lt-s0<(l3t)++(^t-80<(l3[)-]dA  . 


Equation  (3*26)  states  that  the  lnction  in  the  x-direeticn  is 
a small  oscillation  relative  to  a motion  with  uniform  speed 
38  cons' t,  Equation  (3,27)  is  an  expression  of  Archimedes' 
law:  the  rest  position  of  equilibrium  must  be  such  that  the  weight 
of  the  water  displaced  by  the  ship  just  equals  the  weight  of  the 
ship.  The  center  cf  buoyancy  of  our  ship  is  in  the  plane  of  sym- 
metry, and  equation  (3,28)  is  an  expression  of  the  second  law  cf 
equilibrium  of  a floating  body;  namely  that  the  center  of  buoyancy 
for  the  equilibrium  position  is  on  the  some  vertical  line,  the 
y'-axis,  as  the  center  of  gravity  cf  the  ship. 

The  function  ^ must  satisfy 

* ‘‘iTV  * *l„  “ 0 

in  the  domain  D f A where  D is  the  half  space  y < 0,  and  A is  the 
plane  disk  defined  by  the  projection  of  the  submerged  hull  on  the 
xy-plano  when  the  ship  is  in  the  equilibrium  position#  We  assume 
that  A Intersects  the  xz-plone.  The  boundary  conditions  at  each 
side  of  A are 


+ 
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3it 


(3.36), 


u. 


lz 


»o(V*21)+(tol+#21)x  • ®U(y^) 

"flo(hx+®2i^  + <wi+*2i)x  - «n(y-y£> 


The  boundary  condition  at  y = is  found  by  sliminating  n^  from 
(3*10)  and  (3.11).  Since  u>Q  « 0 these  equations  are 


S’h  + »<Ax  - ''it  * 0 
^ly  ’ "o^Lx  + nxt  ■ 0 

* 

and.  they  yield 

(3.37)  4*rxx  * 2»o*lxt  + «ly  *ltt  * 0 

for  y * 0.  The  boundary  conditions  (3*36)  and  (3»37)  show  that 
depends  on  cj^tt),  ©^(t)  and  T*1©  potential  problem 

can  theoretically  be  solved  for  in  the  form 

^1  85  ^ifx*y***ti  W1  ( ^ ) » ®n ( ^ ) » ®21  ( ^ ) ) 


without  using  (3*29),  (3*3°)#  (3»32).  The  significance  of  this 
has  already  been  discussed  in  Seo*  1 in  relation  to  equation  (1.14). 
Hie  general  procedure  to  be  followed  in  solving  all  problems  was 
all  discussed  there* 
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The  general  potential  problem  defined  above  will  be  the  sub- 
ject of  a separate  study.  The  remainder  of  this  paper  is  concerned 
with  the  special  case  cf  a ship  which  moves  along  a straight  course 
into  waves  whose  crests  are  at  right  angles  to  the  course.  For 
this  case  there  are  surging,  heaving  and  pitching  motions,  but 
©!  » 0,  ©2  * 0*  0 and  the  potential  function  / is  an  even 

function  of  z.  Under  these  conditions  the  equations  of  motion  are 
mere  simple.  They  are 

(3.38)  M1S1  = 2p/hx(/lt-s0/lx)dA  + T 


(3.39) 

(3.14-0) 


lyi  = -Spgy^hdx-Spge ^/xhdx+Sp/h  (/lt-sQ/lx)d; 
L L A J 

I31*®31  * -2fS®31/(y-yphdA-2pgy1/xhdx 

A L 

-2pg©31/x2hdx  +iT 


+2p/[x hy-(y-y«)hx](^lt-so/lx)dA. 


It  will  be  shown  in  the  next  section  that  an  explicit  integral 
representation  can  be  found  for  tho  corresponding  potential  function 
and  that  this  leads  tc  integral  representations  for  s^,y^  and 


i 

- h- 
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4*  Method  of  solution  of  the  problem  of  pitching  and  heaving  of  a 
ship  In  a flea-way  having  normal  Incidence. 

In  this  section  we  derive  a method  of  solution  of  the  problem 
of  calculating  the  pitching,  surging,  and  heaving  motions  in  a sea- 
way consisting  of  a train  of  waves  moving  at  right  angles  to  the 
course  of  the  ship,  which  is  assumed  to  be  a straight  line  (i»e. 
u>  a 0).  The  propeller  thrust  is  assumed  to  be  a constant  vector. 

The  harmonic  function  ^ and  the  surface  elevation  ni 
therefore  satisfy  the  following  free . surf ace  conditions  (cf.  (3,10) 
and  (3.11),  with  d>  * 0): 

6 'll  + ao*lx  - +lt  = 0 

(4*1 ) at  y « 0 , 

*ly  * V‘lx  * nlt  * 0 

The  kinematic  condition  arising  from  the  hull  of  the  ship  is 
(cf.  (3.14)  with  e21  * eu  * Wj  » 0): 

(4.2)  *1,  * >0hx  • 

Before  writing  down  other  conditions,  including  conditions 
at  oo,  we  express  ^ as  a sum  of  two  harmonic  functions,  a3  follows: 

(4.3)  (x, y, 8 * t ) * %yU,7,z)  ♦ 


- .as  - 


Here  is  a harmonic  function  independent  of  t vfaioh  is  also  an 
even  function  of  z*  We  now  suppose  that  the  motion  of  the  ship 
is  a steady  simple  harmonic  motion  in  the  time  when  observed  from 
the  moving  coordinate  system  o-x,y,z*  (Presumably  such  a state 
would  result  after  a long  time  upon  starting  from  rest  under  a 
constant  propeller  thrust*)  Consequently  we  interpret  ^3.(x'y»z) 
as  the  disturbance  caused  by  the  ship,  which  therefore  dies  out  at 
oo ; while  7^(x,y,zjt)  represents  a train  of  simple  harmonic  piano 
waves  covering  the  whole  surface  of  the  water.  Thus  is  given, 
with  respect  to  the  fixed  coordinate  system  0-X,Y, Z,  by  the  familiar 
formula 

-2 

2—Y  JcL 

X,  a C e s sin  <dt  + ~ X + p). 


with  6 the  frequency  of  the  waves*  In  the  o-x,y, z system  we  have, 
therefore : 


£v 

X1  ( x, y , z j t ) * C e 8 


sin 


x + (6  + 


+ p] 


We  observe  that  the  frequency,  relative,  to  the  ship,  is  increased 
above  the  value  6 if  sQ  is  positive  **  i*o.  if  the  ship  is  heading 
into  the  waves  - and  this  is,  of  course,  to  be  expected.  With  this 
choice  of  Xj..  It  is  easy  to  verify  that  X0  satisfies  the  following 
conditions: 
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(4.5) 


B^X  + gX 

o*oxx  a*"oy 


* 0 


at  y * 0 , 


obtained  after  eliminating  ^ from  and 


(4.6) 


<oz  * 8ohx  on  A' 


with  A , as  above,  the  projection  of  the  ship's  hull  (for  z > 0) 
on  its  vertical  mid-section.  In  addition,  we  require  that  0 

at  oo  • 

It  should  be  remarked  at  this  point  that  the  classical 
problem  concerning  the  waves  created  by  the  hull  of  a ship,  first 
treated  by  Mchell  [8],  Havelock  [2],  and  many  others,  is  exactly 
the  problem  of  determining  XQ  from  the  conditions  (4*5)  and  (4«6). 
Afterwards,  the  insertion  of  3 XQ  in  (3.33 )»  with  * 0, 

■ 0,  leads  to  the  formula  for  the  wave  resistance  of  the  ship 
i.e.  the  propeller  thrust  T is  determined.  Since  y^  and  0^  are 
independent  of  the  time  in  this  case,  one  sees  that  the  other 
dynamical  equations,  (3*39)  and  (3*4^)»  yield  tho  displacement  of 
the  c.g.  relative  to  the  rest  position  of  equilibrium  (the  so-called 
heave),  and  the  longitudinal  tilt  angle  (called  the  pitching  angle). 
However,  in  the  literature  cited,  the  latter  two  quantities  seem 
to  be  taken  as  zero,  which  implies  that  appropriate  construing 
forces  would  bo  needed  to  hold  the  ship  in  such  a position  relative 
to  the  water.  However,  the  main  quantity  of  interest  is  the  wave 
resistance,  and  it  is  not  affected  (in  the  first  order  theory,  at 
least)  by  the  heave  and  pitch. 


>5>6 


We  proceed  to  the  determination  of  ) { t using  a method  differ- 
ent from  tho  classical  method  and  following,  rather,  a course  which 
it  is  hoped  can  be  generalized  in  such  a way  as  to  yield  solutions 
in  more  difficult  cases* 

Suppose  that  we  know  the  Groen*s  function  G (|,*i,4;x,y,z) 
such  that  G is  a harmonic  function  for  0,  % > 0 except  at 

(x,y,z)  where  it  has  the  singularity  l/r;  and  G**  s tisfies  the 
boundary  conditions 

(4.7)  G^  + kG^  * 0 on  tv  * 0 

G^  = 0 on  K = 0 

where  k = g/s^  • Let  2:  denote  the  half -plane  ^ * 0,  4 > 0; 

and  let  $2  denote  the  half -plane  4*0,  \<  0.  Green*  s formula  shows 

that 


knXo  - - IJX ♦ //  *0-  »*<»««  -ff  Xo40*d|dn  . 
2:  2:  Q 


Then,  since 


- //  ♦//  V*dW  ' ^ ^ lX-°U  ' 

* If  (Xo^  * 


I 
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we  have  an  explicit  representation  of  the  solution  in  the  form 

*°(W)  a^//VVvj  » or 

SI 

(4.8)  ■X0(x,y,S)  - S jj 2.  jj  h{U,>,)0*(?,,,  c ix, y,i)d«<J^ 

A 

upon  using  (k»6). 

In  order  to  determine  Q*  consider  the  Green* s function 
y,z)  for  the  half  space  n<  0 Which  satisfies 

+ kGtt  * o 

on  >1-0.  This  function  can  be  written  as 

G = i • I 

rl  r2  + 8 

where 


V (?  -x ) 2+  (^-7 ) *+«-z ) 2 


1 


— 1 

V (t-x)*+(>i+y )*♦(£„*  )2 


>$8 


and  g is  a potential  function  in  ^ < 0 which  satisfies 


su  tk*n.  = 2k*y 


Vlt-xJ^+y^+di-z)2 


on  \\a  0*  The  well-known  formula 


2k  =§- 

hy 


V Cg-x)^+y2+U-*)* 


2k  J pepyJ0[p^-x)2+(£-z)2]dp, 


in  which  the  Bessel  function  J can  be  expressed  as 

o 

vl2 

JQ  [ p V(  £ -x ) *>  ( £-z ) ^ ] § J cos[pU-x)cos  ®]ccs[p(5-z )sin©]d®, 

o 


allows  us  tc  write 

oo  n/2 

+ kgn  * ¥f  / pepy  cos[p(4 -x)ccs®]cos[p(£-z)aln®]d®dp 

* o o 


for  ^*0  and  y •<  0.  It  is  now  easy  to  see  that 

oo  rr/2 

gff  + kg>.  * ¥//  peP^y+^0os[p(£-x)cos®[p(S-z)sine]dedp 

o c 


is  a potential  function  In  ^ < 0 whioh  satisfies  the  boundary 
condition*  An  interchange  of  the  order  of  integration  gives 


, \ 


-2s*  . * ' 
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ff/2  oo 

/ detLj  P «o.tp(«-*).n.  #1.P[  (yn)*i«-x)c°s#] 


where  ^denotes  the  real  part.  If  we  think  of  p as  a complex 
variable,  the  path  from  0 to  oo  in  the  last  result  can  be  replaced 
by  an  equivalent  path  L: 
rr/2 

i,v  r r p[  (y+«i)+i(5-x)cos  e] 

g^+kg^  = w J d®”t  J p 003  Ip(^”z)®in®3e  dp  . 


Since  the  right  hand  side  of  this  differential  equation  for  g is 
expressed  as  a superposition  of  exponentials  in  % and  n,  and 
since  some  freedom  is  allowed  in  the  choice  of  L,  it  is  evident 
that 

p[  (y+>*)+i(£-x)cos  e] 

g * k \ [ d©  fix  J P o?8 [Baffin  SLe  dP 

™ ’ T kp-pZC032  9 


provided  L can  be  properly  chosen.  The  path  L,  which  will  be 
fixed  by  a condition  given  below,  must,  of  course,  avoid  the  polo 
at  p * k/cosZ©. 

It  can  new  be  seen  that  the  function  0 (£/*7,£fx, y,z)  * 
G($#»l»Six#y,z)  + y»-2)  satisfies  all  of  the  conditions 

imposed  on  the  Green*  s function  employed  in  (lj..8):  the  sum  on  the 
right  has  the  proper  singularity  in  &> 0,  it  satisfies  the 

boundary  condition  (I4..7)  and 


i 
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^U,n»Kix,y,z)  + y#-z) 


is  zero  at  £ ■ 0*  Therefore 


0 ]«-0  * 2 


[v  (*-x)2+(n- 


+(n-y)2+*^  V ($-x)^+(v*+ y)^+z2 


tt/2 


+ f£  /d®^/  cos(pzsinO) 


p[  (y+n)+i(?-x)cos  ©] 
e dp 


k-p  cos  9 


The  substitution  of  this  in  (4.8)  gives  finally 


Xc(x,y,?)  =*  -£f/h  — A_—  - '-V-'-— 

A (VU-x)2+U-y)2+z2  V(5-x)2+(n+y) 


?+z2 


d?d^. 


A C 


tt/2  r p[  (y+V{)+i(g  -x)cos©]  % 

0jJbr((AAj  aelt>jccafPz  3lr‘|j°  — 3e_ /dfd?. 


k-p  ccs  © 


A condition  inposed  onX0(x,y,z)  is  that  X.  (x,y,z)  *■  o as  x +x  • 
This  condition  is  satisfied  if  we  take  L to  be  the  path  shewn  in 
Pig*  4*1 • 


fc 


Fig*  4*1. 

The  Path  L in  the  p-plano 


The  function  is  given  by 


K = *1 + '4  = c 6 g sinnr + 


2-X  2 s,C2  ~1 

8 sin  ^ + (64X.)t+pUt0 


end  therefore 


(4.*9)  /lt  - ®0^ix  * c*e  cos 


»,*2  1 

♦ <«+— T— )t+p  - . V 

o J - c 


If  this  is  substituted  in  the  equation  fer  the  surge  we  have 
• r 2 s 62  "j 

M^s-^  * 2pC6  JJ  hxo  6 c08  £-2  (6+  ) t+p  dxdy 

- 2ps„  //  h X dx  dy  + T • 


The  last  equation  shows  that  in  order  to  keep  s1  bounded  for  all  t 
we  must  take 


(4.10) 


T - 2pa0  JJ  hx  Hox  dxdy 


where 


XO3C(x.y,0)-.^jr/h  (8^) 


LjfcsL 

4«-x)2+(tl-; 


y)2]3/a  [(f,x)c+(Yi+y)c] 


2s, 


// 


f”/2 

V*,Ylj/  c 


pC  (y+nJ+K'C-xJcos©] 


/ d9%  ] hs,~ssti  e 5 2E. 

L S“so  P c°s  « 


dfd**. 


Equation  (4*1°)  gives  the  thrust  necessary  to  maintain  the  speed  sQ, 
or  inversely  it  gives.:the  speed  so  which  corresponds  to  a given 
thrust*  The  integral  in  (4*10)  is  called  the  wave  resistance 
integral*  As  one  sees,  it  does  not  depend  on  the  seaway*  The 
integral  can  be  expressed  in  a more  simple  form  as  follows* 

The  function  ^ox(x,y, o)  is  a sum  of  integrals  of  the  type 

//^(5,n)  dfdn* 

A 

If  an  integral  of  this  type  is  substituted  in  the  wave  resistance 

a* 

integral  we  have 


*$*v.r**»* 


JJ  SI  hx(x,y)h-(l,n)f(«^ix»y)  dfdttfxdy 

A A 

say*  This  Is  the  sane  as 

JJ  ff  h^(«,n)hx(x,y)f(x,y|t#n)  dxdydfcdV 

A A 


and  we  see  that  I = 0 if 


f(£#^x,y)  * - fU,y }£,yO  * 


Therefore 


T **  JJ  JJ  hx(x,y)h£(£,*i)  tl  dfcdrjdxdy 


A A 


where 


p(y+n) 

cos©  e co 


g - sQ  p cos  © 


IV 


Since  & / is  zero  except  for  tho  residue  from  the  integration  along 


the  semi-circular  path  about 


8.008  0 cos  © 


we  find  from  the  evaluation  of  this  residue  that 

, k(y+n)  sec2© 

fj  * 2|L-  j secJ9  e cos  [k(g-x)  oos  9]  dfi  • 

so  o 

Now  if  we  define 

2 

P ( 0 ) * jj  hx(x,y)  e1^  sec  ® cos  (kx  sec  9)  dxdy 
A 

2 

Q(9)  * Jf  hx(x,y)  o*7  sec  6 sin  (kx  sec  9)  dxdy 
A 

we  can  write 

2 */2 

T » Jttg-  J (p2  + Q2)  see3  9 d9  # 
nso  o 

This  is  the  familiar  formula  of  Ilichell  for  the  wave  resistance. 
The  surge  is  givun  by 


Hereafter  we  will  suppose  for  simplicity  that  there  is  no 

coupling  between  (3.39)  and  ( 3 , sothat  / xhdx  * 0,  The 

L 

substitution  of  (I4..9)  in  (3*39)  there  gives  the  following  equation 


for  the  heave# 


Mjy-j+fapg  f hdx]y1 
L 


6y 
r r ■"■j** 

,2  S-62 

//h  e 6 cos 

V+(s+i-,ttp 

¥ 

A 

dxdy 


-2P  so  //  V*  ixA* 

A 


The  time  independent  part  of  y^,  the  heave  component  of  the  trim, 
we  denote  by  y^j  it  is  given  by 

(4*  11)  (g  J h dx)y*  * -sQ  //Vox  dxdy. 

L A 

y^  is  the  vertical  displacement  of  the  center  of  gravity  of  a ship 
moving  in  calm  water  from  its  rest  position.  The  integral  on  the 
right  hand  side  of  (4.11)  is  even  more  difficult  to  evaluate  than 
the  wave  resistance  integral.  As  far*  as  the  authors  are  aware, 
the  Integral  has  not  appeared  in  the  literature. 

The  response  of  y1  to  the  sea  is  given  by 


p r 0 

fA2v  f2  1 

2pCo 

• 

i 

[ j hye  6 cos 

i 

£-5  + (o+s  ~)  t+p 
g v o g'  * 

“ m 

dxdy 

f °r6  ? 

2pg  J hdx  - )* 

L 6 


For  the  case  under  consideration,  the  theory  predicts  that  resonance 
in  the  heave  occurs  when 


The  equation  for  the  pitching  angle  ia 


I31°31+2ps 


[ / ly- 


yJ)hdA  + / x2hdx  | S 


31 


2pC<5"  jj  [xhy-(tr-y^)hx]cos|~^  +($*•*-£ — )t+p|  dxdy 


+ iT  - 2ps0  /[xh  -(y-y*)hx]  Xox  dA  . 


The  time  independent  part  of  9^,  9^,  is  given  by 
2pg  £ J (y-y* )hdA  + / x2h  dxj  9^1 

**  i^T  - 2ps0  / [xhy-(y-y»  )hx]  Xox  dA  . 

a ( i-y«)T  - 2ps0  y [xhy-yhxJXox  dA# 

A 

& 

The  angle  is  called  the  angle  of  trim;  it  is  the  angular  dis 
placement  of  a ship  which  moves  with  the  speed  s in  calm  water. 
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